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Abstract: We study the implications and limitations of galaxy cluster surveys for con- 
straining models of particle physics and gravity beyond the Standard Model. Flux limited 
cluster counts probe the history of large scale structure formation in the universe, and as such 
provide useful constraints on cosmological parameters. As a result of uncertainties in some 
aspects of cluster dynamics, cluster surveys are currently more useful for analyzing physics 
that would affect the formation of structure than physics that would modify the appearance 
of clusters. As an example we consider the ACDM cosmology and dimming mechanisms, such 
as photon-axion mixing. 
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1. Introduction 

In this era of precision cosmology, a wide variety of cosmological and astrophysical observa- 
tions are providing strong constraints on the composition of our universe. Among these are 
studies of the cosmic microwave background (CMB) Q, large scale structure luminosity- 
redshift curves of Type la supernovae ||, galaxy rotation curves and light element abun- 
dances Q . A relatively consistent picture of the universe has emerged in which the current 
universe is flat (f2 = 1), contains about 20% of its energy density in nearly pressureless cold 
dark matter, about 76% in dark energy (Q\ = .76), and the remainder in ordinary matter 
described by the Standard Model of Particle Physics (SM) jj]]. (We take Q m to be the sum 
of the cold dark matter and Standard Model matter, including neutrinos, so O m = 0.24 
by the above estimates.) The flatness of the universe and the spectrum of initial density 
perturbations is explained by the paradigm of inflation. On the other hand, dark matter 
and dark energy provide a challenge for particle physics. The influence of dark matter on 
galaxy rotation curves, the CMB, and most directly in the observed separation of dark matter 
and baryonic matter in the "Bullet cluster" ||, provides conclusive evidence that there are 
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new types of particles which have not been observed in particle physics experiments and are 
not described by the SM; and it may be argued that the observation of dark energy in the 
expansion history of the universe hints at new gravitational physics. 

The incredible precision of lunar ranging measurements produce some of the strongest 
constraints on new gravitational physics but only on local phenomena that would be 
occuring here and now. The overall expansion history of the universe constrains the influence 
of new physics on the largest of scales, and indeed the luminosity-redshift curves of Type la 
supernovae have provided the most direct evidence for dark energy. The formation of structure 
in the universe is also highly dependent on gravitational and particle interactions, and since 
structure has had a relatively long time to form, galaxy and cluster surveys provide another 
useful probe of the amount and features of dark matter and dark energy, as well as other new 
physics. The purpose of this paper is to examine the importance of galaxy cluster surveys in 
testing of new ideas in gravitational and particle physics. (See also || p.) It is certainly not 
a new idea to use structure formation to constrain cosmological models. Indeed, the Press- 
Schechter formalism for predicting counts of virialized objects is more than 30 years old fi~0[| . 
Clusters are the largest virialized objects in the universe, and as such provide a useful probe 
of structure formation. Collisions of hydrogen atoms in the intracluster gas produce X-rays, 
and track the gravitational potential well in a cluster [11|. As a result, X-ray surveys have 
provided reliable and complete surveys of galaxy clusters in various regions of the sky. Serious 



studies of the properties of X-ray clusters for this purpose began in the 1980's |12|, [13|, [b|, 15| . 
It was suggested by some groups that cluster surveys were in conflict with the Concordance 
Model (Q m = 0.3, a 8 ~ 0.9) gf| 0, H, H. It now seems that the HIFLUGCS cluster survey 
is in agreement with the most recent Cosmic Microwave Background data (fi m = 0.234 and 
cr 8 = 0.76) 0, 1^. Our analysis with the ROSAT 400 Square De gree data set is also in 
agreement with CMB data if a large scatter is assumed in the relation between cluster mass 
and temperature. Otherwise, our analysis prefers cosmological parameters closer to the old 
Concordance Model. 

The particle physics community has not yet embraced cluster technology for the purpose 
of testing physics beyond the standard model. In large part this is because of limited statis- 
tics and uncertainties in the theoretical models of structure formation and cluster dynamics. 
However, while supernovae are sensitive to the geometry (and opacity) of the universe, while 
structure growth is sensitive to its clustering properties, these are truly complementary ap- 
proaches, as argued by Wang et. al [21]. However, since cluster surveys can provide constraints 
on new physics complementary to other cosmological constraints, they deserve to be in the 
arsenal of the particle physics trade. A purpose of this paper is to review and introduce 
much of the technology involved to the particle physics community. As an example of the 
application of cluster surveys to particle physics and its limitations, we study the significance 
of current and future surveys for constraining dimming mechanisms such as the photon-axion 
oscillation model of Csaki, Kaloper and Terning [22|. To motivate consideration of dimming 
mechanisms, we note that while there are numerous models of particle physics beyond the SM 
which provide dark matter candidates, the nature of the dark energy is more of a mystery. 
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Constraints on the dark energy equation of state from WMAP and the Supernova Legacy 
Survey (SNLS |23| ) suggest that, assuming a flat universe, w = — 0.97±0.07 Q, where p = wp 
is the linearized equation of state relating the pressure of the dark energy fluid p to its energy 
density p. The value w = —1 describes the vacuum energy, or cosmological constant. How- 
ever, naive particle physics estimates of the vacuum energy are dozens of orders of magnitude 
too large, so it is well motivated to consider alternative models. 

If a new pseudoscalar particle existed with a certain range of mass and axion-type coupling 
to the electromagnetic field, then distant objects would appear dimmer than expected because 
a fraction of the light emitted by the stars in a galaxy would have been converted to axions 
while traversing the intergalactic magnetic field |2^] . It would be necessary to reevaluate 
the evidence for acceleration of the universe if the dimming of distant supernovae could be 
explained without a cosmic acceleration. At the time when the photon-axion oscillation model 
was proposed, a universe without acceleration could not be ruled out if one allowed for such 
a dimming mechanism. Since that time, new data has provided stronger constraints on the 
equation of state parameter of the dark energy, and a model without acceleration is currently 
disfavored Q. However, photon-axion oscillations (or any other viable dimming mechanism) 
could still exist, and would lead to an apparent decrease in the dark energy equation of state 
parameter, a possibility which remains open [24]. 

A dimming mechanism would also affect flux limited galaxy cluster surveys. Some distant 
galaxy clusters which would have otherwise been bright enough to be detected in a flux limited 
telescope, may become too dim to be detected as a result of photon loss, but it is not a priori 
obvious what the implications for cosmological analyses would be. Although the appearance 
of clusters would be affected by such dimming, such effects can be absorbed into the measured 
evolution of the luminosity-temperature relation. At any rate, clusters provide an independent 
test of the nature of dark energy which is complementary to supernovae, and thus potentially 
constraining of models such as dimming mechanisms. 

In the following, we will attempt to provide a thorough review of how cosmology re- 
lates to theories of structure formation and our observations. Our analysis relies on several 
assumptions regarding cluster luminosities and their evolution. A better theoretical under- 
standing of the evolution of cluster properties is desirable (however, see Ref. |2^]). On the 



other hand, since the apparent evolution of cluster luminosities has been measured |26|, 27], 
cluster counts provide more direct constraints on new physics that would affect the formation 
of clusters rather than their appearance. We will use the 400 Square Degree ROSAT survey 
]pg[] (hereafter referred to as 400d) as our primary data set. We also use the 400d survey to 
constrain the standard ACDM cosmology, which does not require a theoretical understanding 
of the mechanisms of luminosity evolution. 

In Section Q we review the statistical models of structure formation based on the Press- 
Schechter formalism |10[| . In Section |3] we analyze the possibility of photon-axion oscillations 
in light of current galaxy cluster surveys. Interestingly, while supernovae surveys can be dra- 
matically affected by dimming, because the redshift evolution of luminosity and temperature 
is measured, the studies of cluster count evolution are remarkably insensitive to it (although 
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the total counts, themselves, are). In Section || we present our statistical analysis of cluster 
constraints for standard cosmology, and thus demonstrate the techniques which are simply 
applied to other theories of modified dark matter or dark energy. In Section || we examine 
the significance of future cluster surveys for probing new gravitational and particle physics. 
We conclude in Section || 



2. Analytical Models of Structure Formation 



Although the state of the art in structure formation involves elaborate n-body simulations, 
much can be understood within simple, analytical models. In this section we summarize the 
basic theory behind structure formation in the universe and how the theory is compared with 
galaxy cluster surveys. The review is simplified, and does not contain new results, but we 
hope it contains enough of the basic ideas so that particle physicists can easily apply the 
formalism to constrain new physics. There are a number of excellent reviews on structure 
formation in the literature that are substantially more comprehensive than this one, such 
as Refs. 2S]. Techniques for comparing the models to X-ray cluster data are somewhat 
scattered in the literature, though recent cluster surveys provide useful background with their 
catalogues, as in Refs. [17, 2£, j3Cfl . Our goal is to simplify the discussion to its bare essentials 
without forfeiting too much of the underlying physics, making use of the fact that others 
have performed the complicated simulations necessary to test both the phenomenological 
models of hierarchical collapse, and hydrodynamic scaling relations between cluster mass and 
observational quantities like cluster temperature. Simulations suggest that the simplified 
models of structure formation and X-ray cluster dynamics are accurate enough to constrain 
new physics by building the new physics on top of these models. While simulations are not in 
perfect agreement with these models agreement is good enough that these models serve 
as a useful tool in studying the evolution of structure. 



2.1 The Press-Schechter Formalism 

The CMB provides strong evidence that the universe was homogenous to a part in 10 5 at the 
time that atoms formed during recombination. However, as the universe expanded structure 
formed due to the gravitational collapse of these small fluctuations into progressively larger 
objects. The precise way in which structure formation occurs is sensitive to the composition 
of the universe. Smooth, unclustering dark energy, for example, leads to a faster expansion 
of the universe and hinders formation of structure on large scales. Since the evolution of 
structure depends on the composition of the universe, comparison of models to observations 
provides an important probe of cosmology. The Press-Schechter (PS) formalism provides 
a simple model for translating cosmology into number counts for structures on arbitrary 
length scales, as a function of mass and redshift. Here we summarize only the main results 
of this formalism, but there are many lengthier discussions in the literature justifying this 
approach and deriving the relevant results quoted below (e.g. Refs. [f29f). 
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The spherical collapse model of Press and Schechter imagines that initially overdense 
regions of the universe collapse with spherical symmetry, in an otherwise homogeneous uni- 
verse. As the universe expands, only objects with a density above a critical value will have 
collapsed and virialized by any given time. For the remainder of this discussion we assume a 
flat universe with Q m + £l\ = 1. In terms of the cosmological parameter ttt(z), where 

n f (z) = (Q m (1 + zf)/(Q m (l + zf + (1 - J2 m )), (2.1) 

the critical overdensity at the time of virialization in the linearized spherical collapse model 
is given by [32]: 

51 * { 2Q j (1 + 0.0123 log 10 O,) . (2.2) 

In the Einstein-de Sitter universe with Q m = 1, 6g C ~ 1.69. In order to compare with the 
density field observed today, we need to account for the evolution of the universe. Since 
virialized objects are nonlinear fluctuations of the background density field, it would seem 
difficult a priori to describe analytically the evolution of the statistical distibution of density 
perturbations. It was the observation of Press and Schechter that since the original spectrum 
of density perturbations was approximately Gaussian, and because the precise nonlinear evo- 
lution of those density perturbations is unlikely to significantly modify the mass contained 
in collapsed objects, a linearized approach can be justified for modeling the distribution of 
massive collapsed objects. 

In the linear model, density perturbations grow proportional to the growth factor D(z), 
and thus the linearized overdensity of an object that virialized at a redshift z has grown by 

5 c {z) = 5l c D{0)D{z)-\ (2.3) 

where the linear growth factor is defined as [33, 34 1, 



(1 + z') 

D(:) — 2"){'1 I hHqH(z) I ±—J-dz', (2.4) 



and H{z) is the Hubble parameter at redshift z, 



H(z) = H oy /n m (l + zf + n A . (2.5) 

Below we will use H = lOO/i km/s/Mpc with h = 0.73 g]]. 

If one assumes a Gaussian distribution, the probability of a given collapsed object of mass 
M having an overdensity in the linearized model larger than S c (z) today is 

1 f°° 

p(S c (z),M) = -^ exp(-5 2 /2a 2 )d5. (2.6) 

V27R7 J8 c (z) 

By differentiating p(5 c , M) with respect to M and dividing by the volume (M/p) one gets the 
number density of objects with mass between M and M + dM. The present day mean matter 



density of the universe is p = 2.775 x 10 n tt m h 2 M Mpc-^im. In this simplified picture of 



- 5 - 



structure formation, small objects become bound first, followed by larger structures. Galaxy 
clusters are the largest virialized objects in the universe in the current epoch, which makes 
them especially suitable for study by this approach. 

The variance in the distribution of density fluctuations in the universe, a{M) 2 = ((5M/M) 2 ) 
is typically normalized to spheres of radius 8 ft. Mpc: 

, ,2 2 f^k 2+n T(k) 2 \W(kR(M))\ 2 dk , 

a (M) = al — ^ K - L ^ V " ] 9 — (2.7) 

j™k 2 + n T{k) 2 \W{k%h- l M V c)\ 2 dk V ' 

where a% is fit to cosmological data. The function W(x) is a top-hat filter, 
W(x) = 2 (sinx — xcosx) /x 3 , so that a(M) is the variance of the mass distribution in 
spherical volumes of radius R{M) = [3M/ (47rpo)] 1 ^ 3 - The initial power spectrum is usu- 
ally assumed to take the scale-free Harrison-Zel'dovich form P(k) oc k n , with n = 1. As the 
universe expanded density perturbations grew at different rates through the radiation domi- 
nated era to the present era. Likewise different size fluctuations crossed the particle horizon 
at different times. The power spectrum therefore varies from the Harrison-Zel'dovich spec- 
trum in a wavelength dependent way. Assuming linearity, i.e. an absence of mixing of modes 
with different wavenumbers, the power spectrum evolved from early times may be written 
as P(k) = T(k) 2 k, 1 where the transfer function, T(k), takes the primordial power spectrum 
to the present day. As k — > 0, T(k) — > 1 because large enough wavelength fluctuations 
have not crossed the particle horizon and therefore keep their primordial spectrum. This also 
means that the initial time can be taken to be any time before which fluctuations of interest 
would have crossed the horizon. The form of the transfer function is found by analyzing 
numerically physical processes that would modify the power spectrum over time. For cold 
dark matter, Bardeen, et al. [p8| found, 



T(k) = In (1 + 2.34?) / (2.34?) 

x [l + 3.89<7 + (16.1?) 2 + (5.46g) 3 + (6.71g) 4 ]~ 1/4 (2.8) 

where q(k) = kj (p. m h 2 Mpc -1 ) in the absence of baryonic matter. To account for baryon 
density oscillations a "shape parameter" V is introduced, simply replacing q(k) by 



j(k) = -, (2. 

y ' r/iMpc- 1 v 



where [40] 



r = Sl m h exp 



n b (l + V2h/n m )] , (2.10) 



and f^b is the ratio of the baryon density to critical density. 

While the Press-Schechter formalism is remarkably successful in its comparison to numer- 
ical simulations (|l0|, ^2j), it has proven to be most powerful as a basis for a phenomeno- 
logical approach to modeling galaxy cluster counts. One extension to the formalism takes 



1 We absorb the superhorizon evolution into P{k); see Ref. feg for a discussion. 
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into account non-sphericity of collapsing objects. Sheth, Mo and Tormen (SMT) developed 



a modified PS procedure [41] which, allowing for ellipsoidal collapse, introduces new model 
parameters which are fit to N-body simulations. In the SMT model, the mass function is 
given by 

dN [2a~ p dv ( 1 \ / av 2 \ , 

dM = yV C MdM{ 1 + Ja^Yp) eXP ) ' (211) 

where v = 5 c (z)/a(M), and the best fit for the parameters a, c and p assuming a standard 
ACDM cosmology are a = 0.707, c = 0.3222, p = 0.3 |4l|]. (By comparison, in the PS model, 
a = 1, c = 0.5 and p = 0.) 

2.2 Relating Measured Flux to Cluster Luminosity and Mass 

The Press-Schechter formalism and its extensions reviewed above predict the statistical distri- 
bution of massive collapsed objects in the universe as a function of their masses and redshifts. 
On the other hand, telescopes do not directly measure cluster masses, but rather the flux and 
perhaps the spectrum of light emitted by those clusters in some frequency band as observed 
on or near Earth. In order to relate the mass function ( 2.11| ) to observational quantities, it 



is necessary to understand the relationships between the mass of a cluster and observational 
data. In this section we describe how properties of X-ray clusters are related to one another, 
and how those properties are then compared with observations. 

On average, hydro dynamical models which yield simple scaling relations between cluster 
mass and cluster temperature (the M — T relation) have been proven reasonably successful 
in comparison with numerical simulations |35|, [43|]. On the other hand, the relation between 
the temperature and X-ray luminosity (the L — T relation) of clusters is sensitive to more 
complicated physics such as cooling mechanisms and the density profile of the intracluster 
gas, and is fit by cluster data. Furthermore, it is now commonly accepted that the L — T 
relation has evolved as the composition of radiating cluster gases has evolved [gfj] . 

There are at least two sensible notions of cluster temperature, so it is important to be 
precise in terminology. From here on when we refer to a cluster's temperature, T, we will 
mean the temperature of the baryonic gas in the cluster, as is directly measured from the 
spectrum of light emitted by the gas in the cluster. We model the cluster gas as isothermal, 
which may not be that good an approximation for actual clusters [E4|, although predictions 



for number counts are not that sensitive to this assumption [35]. Another notion of cluster 
temperature is determined by the velocity dispersion of the dark matter particles, a 2 = (v 2 ), 
where the velocity v is measured in the rest frame of the cluster and the brackets denote the 
statistical average over dark matter particles. If typical dark matter particles have a mass 
rriD, then the quantity To = mDO 2 /ks is a measure of the temperature of the dark matter in 
the cluster, where fee is Boltzmann's constant. Generally To is not directly related to T, as 
the dark matter is not expected to be in equilibrium with the baryonic matter. However, it 
is often assumed that these temperatures are similar, or at least proportional to one another, 
after which a scaling relation between cluster temperature and cluster mass follows. 
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For an isothermal spherical cluster of dark matter, the density p and velocity dispersion 
a scale with distance from the cluster center r as [45]: 



p(r) 



a 



2TrG N r 2 ' 



(2.12) 



with Newton's constant Gn- As mentioned earlier, the assumption of isothermality may 
not accurately describe the density profile of the halo, which is a subject of intense study. 
A phenomenological density profile which fits better simulations is given by the model of 
Navarro, Frenk and White (NFW) (^6|, in which the density profile takes the form, 



p(r) 



2 ' 



(2.13) 



where p s and r s are model parameters. In the current analysis we assume the isothermal 
profile, Eq. ( 2.12| ), for easier comparison to analytic approximations of scaling relations in the 
literature. 

By considering the evolution of spherical density perturbations, one can estimate the 
density of objects which had just virialized at redshift z. The density of virialized objects 
may be written in terms of A(z), the ratio of the cluster density to the critical density, 
Pcrit = 3H 2 /8irG. Assuming unclustering dark energy, and ordinary CDM, a useful analytic 
approximation to A(z) was given in Ref. [fl?]] for flat ACDM cosmologies: 



A (z) = 18vr 2 + 82(0/ (z) - 1) - 39(0/(z) - 1) 



(2.14) 



where fi/(z) is given by Eq. (2.1). A scaling relation between the velocity dispersion and the 
cluster mass is obtained by approximating the mass of a spherical cluster which virialized 
at redshift z to the mass obtained by integrating Eq. (2.12) to a radius such that the mean 



density is given by p cr u A(z), with the result [47], 

M 2/3 jff(z) 2/3 A(z) l/3 > 



2 

(J 



(2.15) 



Here, H(z) is the redshift-dependent Hubble parameter, Eq. (2.5). Assuming the baryonic 
gas in a cluster has temperature proportional to the dark matter velocity dispersion a 2 , it 
follows that the cluster temperature, T, scales with cluster mass, M, and redshift, z, as in 
Eq. O: 

T ~ M 2 l z H{zf z ^{z) 1 ^. (2.16) 
In practice, simulations are used to determine the constant of proportionality T15 defined 



through [17 1, 



- 15 



h \ 2//3 / O m A (z, $7 r , 



0.73 J 



178 Q f (z) 



1/3 



/ M 
[ 10 15 M m 



2/3 



(1 + *) 



(2.17) 
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where Qf(z) is given by Eq. fl2.1|) , and Mq is the solar mass. The normalization factor 178 is 
approximately the overdensity of a just-virialized object (c./. Eq. ( |2.2D in the linear model). 2 

Different simulations determine a variety of values for T15 , which leads to some ambiguity 
as to the most accurate normalization for the M — T relation. Typical values are T15 ~ 4.8 



keV and T15 5.8 keV [47]. We do fits for various values of T15 to gauge the errors associated 
with the uncertainty in the M — T relation. 

With a relation between cluster temperature and cluster mass in hand, the SMT mass 



function, Eq. ( 2.11] ), can be used to determine how many clusters of a given temperature are 



expected per unit volume of the sky as a function of redshift. However, telescopes often have 
poor spectroscopic resolution, so that in many X-ray cluster surveys it is difficult to accurately 
determine cluster temperatures. Furthermore, telescopes are unable to observe arbitrarily dim 
objects, i.e. they are flux limited. Hence, in order to use the Press-Schechter formalism to 
predict observed number counts of galaxy clusters it is still necessary to relate the cluster 
temperature to observed flux. Such a relationship comes in the form of the elusive L — T 
relation f|, M. There are a number of complications in predicting, and in making practical 



use of, the L—T relations which appear in the literature: («) Surveys often quote fluxes in some 
frequency band, not the bolometric (i.e. total) flux. X-ray telescopes are sensitive to light 
with frequencies of fractions of a keV to tens of keV, though not always in precisely the same 
frequency band, (ii) The measured frequency band is specified in the telescope's reference 
frame, so redshifting of the sources affects the fraction of the total luminosity observed in the 
specified frequency band, (m) There is some scatter in the L — T data (from which L — T 
relations are fitted), which is due in part to a complicated cooling process that takes place in 



many clusters in the central region of the cluster gas [14|. As a result, when possible, some 
surveys remove the central cooling regions when inferring X-ray luminosities, and some do not. 
(iv) In addition, there is relatively strong evidence that the L — T relation has evolved over 
time pq] due to changing cluster environments, (v) Furthermore, when inferring luminosities 
of distant objects from measured fluxes a particular cosmology must be assumed, and the 
assumed cosmology may differ from one quoted evolving L — T relation to another. 

In this paper we focus on the recent 400d ROSAT survey [28], so we will make use of 



published L—T data most easily compared to the cluster luminosities as presented by the 400d 
survey. In particular, the 400d survey quotes X-ray fluxes in the 0.5-2 keV band including 



the central cooling regions. We begin with the L — T relation determined by Markevitch 14] 
from 35 local (z < 0.1) clusters. The fitted power law L — T relation takes the form 



- local 



L °°™ 2 - 4 = A ° ( 6keV 1 ' (2 - 18) 



2 Refs. [17, [l8[ define A(z) as the contrast density with respect to the background density at redshift z. As 
in Ref. j47|, we are defining the contrast density with respect to the critical density, p cr it = 3H 2 /8nG. This 
is the origin of the different scaling relations as written in Ref. and in Refs. Jl7| , ^8| . Physically, they are 
equivalent, assuming an isothermal density profile. 
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with A 6 = (1.71 ± 0.21) 10 44 /i _i erg s -1 , and B = 2.02 ± 0.40, where cooling flows were not 
removed when inferring either luminosities or temperatures 14]. 

To study the redshift dependence of the L — T relation, Vikhlinin, et al. |^] measured the 
X-ray temperature and fluxes of 22 clusters at redshifts 0.4 < z < 1.3 and with temperatures 
between 2 and 14 keV. The luminosity L inferred from the flux F depends on the assumed 
cosmology via, 

L = 4ir d L (n m , z) 2 FK(z), (2.19) 

so the observed redshift dependence of the L — T relation depends on the cosmology. The 
K-correction K(z) will be discussed below. The luminosity distance, cIl, is given by, 



d L (£l r 



l' z cdz' 

+ W (2.20) 



The integral is the comoving distance between the source and the telescope. The extra factor 
of (1 + z) 2 in d\ accounts for the decreased energy per photon from redshifting of the source, 
and the decrease in frequency between photon arrival times, as the universe has expanded. 
To correctly interpret the luminosity evolution in different cosmologies, Eq. ( 2.1 8| ) should be 
modified, assuming power law evolution, with the reference cosmology factored out: 

^0.1-2.4 = A 6 (-l-) B (1 + *)« (2.21) 



,6keVy d L (l,z) 2 

Vikhlinin, et al. found that assuming a £l m = 1, S7a = reference cosmology, a = 0.6 ± 
0.3. It is important to stress that a nonvanshing power a does not in itself imply an evolution 
of cluster properties, because a is cosmology dependent. However, assuming a more realistic 
Vt m = 0.3, U\ = 0.7 reference cosmology leads to a still larger power, an ra =o.3 = 1.5 ± 0.3 
Hence, it seems difficult to argue that the inferred luminosity evolution is due to a 
mistaken assumption about the cosmological expansion rate. We also note that other surveys 
find similar results. For example, the XMM Omega project determined a = 0.65 ± 0.21 [17]. 



We need to be able to convert between frequency bands both as a result of the redshift- 
ing of the spectrum, and in order to compare measurements of surveys in different frequency 
bands. From Eq. fl2.21|) , we can infer a similar relation for the luminosity in the 0.5-2 keV band 
(as appropriate for the 400d survey) in the cluster rest frame if we know the X-ray spectrum. 
The difficulty is that spikes in the spectrum from atomic excitations contribute significantly 
to the flux, so some understanding of the components of the cluster gas is necessary to accu- 
rately convert luminosity in one frequency band to luminosity in another frequency band. A 
popular and accurate model is the optically thin plasma model of Mewe, Kaastra, Liedahl, 
and collaborators p9], the so-called Mekal model. The Chandra Interactive Analysis of Ob- 
servations (CIAO) software package [|5(]] contains code for the purpose of converting spectra 
between bands and between reference frames, and includes packaged spectral models. The 
REFLEX cluster survey has also tabulated conversion factors between luminosities in various 
frequency bands [ 30 1 as a function of temperature, for easy comparison of cluster data to 
structure formation models without necessitating installation of the CIAO software. 
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We assume that the luminosity evolution parametrized by the (1 + z) a dependence in 
Eq. ( 2.21| ) is uniform across the spectrum, so that the same power a will describe evolution of 



the 0.1-2.4 keV L — T relation as in the 0.5-2 keV L — T relation. As a test of this assumption 



we studied the z > 0.4 cluster data by Vikhlinin et al. [26], which includes measurements 
of flux in the 0.5-2 keV band and bolometric flux. We checked that Vikhlinin et al.'s fit of 
a ~ 0.6 is valid both with their measured bolometric fluxes and fluxes in the 0.5-2 keV band. 
For our fits we use, 

a = 0.6 (with an £l m = 1 reference cosmology) 
A 6 = 1.06 (for the 0.5 - 2 keV band) (2.22) 
B = 2.02 



as the parameters in Eq. (2.21). 

The final factor required to compare intrinsic luminosity to measured flux is the K- 
correction. The K-correction converts from the luminosity in a specified frequency band in 
the rest frame to luminosity in the same frequency band in the lab frame, as per Eq. ( j2.19| ). 
The K-correction for a source at redshift z in the frequency band (/i, $2) is given by, 

f! 2 dfP T (f) 

K(T, z) = — . ( 2 - 23 ) 

where Pr(f) is the rest frame spectral distribution for an X-ray cluster with temperature T, 
as a function of frequency /. For example, the measured flux F in the 0.5-2 keV band from 
a cluster at redshift z with rest frame luminosity -Lo.5-2 and temperature T is given by, 

p = A).5-2 (2 24) 

ATrd L {n m , z) 2 fl 5 dfp T (f) 

The K-corrections are not strongly temperature dependent except at the low end of typical 
cluster X-ray temperatures, and a simple power law spectrum, 

Mf) ~ r*. (2-25) 



with index n = 0.5, is found to give a reasonable fit in the relevant frequency bands [51]. 
A comparison of the K-corrections from the simple power spectrum and from more precise 
plasma spectra can be found in Ref. pi] ], or from the documentation for the Sherpa module 
of the CIAO software [ |52[] . Since one of our goals is simplicity in comparison of models of 
new physics to cluster data we will assume the simple power spectrum in our fits. One should 
keep in mind, however, that without much more difficulty more accurate K-corrections can 
be obtained using available software. 

Finally, in order to predict the number of observed clusters it is necessary to know the 
probability of the given telescope detecting a cluster with a given flux. The selection function 
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measures this probability, and depends on the particular survey. The selection function is 
often presented as an effective sky coverage area as a function of flux, but can easily be 
converted to a detection probability. The ROSAT 400d survey contained a geometric survey 
area of A geo = 446.3 deg 2 p8fl . The selection probability is obtained from their tabulated 
effective sky coverage, A e ft(f) as a function of flux /, via, 

Psei(f) = (2.26) 

We are now prepared to calculate the number of clusters we expect to see in an area 
of the sky per unit redshift, for a given telescope flux limit. The mass function Eq. ( 2.11| ) 



gives the distribution of clusters as a function of mass and redshift. We convert mass to 
temperature using the M — T relation, Eq. fl2,17| ); and then temperature to measured flux 
in the appropriate frequency band using the L — T relation, Eq. fl2.18|) , K-corrected as in 
Eq. fl2.23| ), with parameters specified by Eq. ( p. 22 ). 



3. Dimming Mechanisms and Cluster Counts 
3.1 CKT photon loss mechanism 

As an example of a dimming mechanism we will study the possibility of photon loss due 
to photon-axion oscillations (the CKT model j22|) as hypothesized by Csaki, Kaloper and 
Terning. We first note that both cluster counts and Type la supernova surveys extend to 
comparable redshifts z > 1. Hence, if the parameters of the CKT model are chosen so as to 
affect the interpretation of the supernova data, as in |22| |, then the same dimming mechanism 
will have an affect on flux limited cluster surveys. 

The CKT model assumes an axion-like interaction between a pseudoscalar field <f){x) and 
the electromagnetic field of the form, 



Hnt 



Ml c 2 



>E B, (3.1) 



where the dimensionful scale Ml governs the strength of the interaction. 

There is significant evidence for an intergalactic magnetic field (IGMF), although little 
is known about its uniformity in magnitude and direction [p3| . The typical magnitude of the 
IGMF is estimated to be 10~ 9 Gauss, and it is typically assumed that the field is coherent 
to about Ldom ~ 1 Mpc. In a background magnetic field the photon-axion interaction, 



Eq. (3.1), gives rise to a mixing between the axion 4> and the electric field. This mixing 
causes oscillations just as mixing between neutrino flavors leads to neutrino oscillations. In 
the case of photon-axion oscillations there is one polarization of the axion and two of the 
photon, so after traversing enough regions of randomly oriented magnetic fields a beam of 
photons will become distributed equally among the three polarizations. Asymptotically, the 
intensity of light received from a distant astrophysical object will be decreased by a factor of 
1/3 (in the limit of infinite horizon size). As discussed in |2^[ , this effect is approximately 



- 12 - 



described by the following expression for the probability of a photon to remain a photon after 
traversing a comoving distance r(z): 



where, 



P^{z)oL\ + \e-rW L *~, (3.2) 



_ 8 h 2 c 6 Ml 

deC ~ ^7 Fd|2 ' ^ ' 



The existence of a dimming mechanism like photon-axion oscillations would modify the in- 
terpretation of cluster counts. In particular, there would be a reduction in the number of 
distant visible clusters in a flux limited observation. Dimming can also mimic luminosity evo- 
lution, although photon-axion oscillations cannot explain the observed luminosity evolution. 



The fact that a > in Eq. ( 2.21 ) indicates that distant clusters appear more luminous than 
nearby clusters, while photon-axion oscillations would have led to the opposite conclusion. If 
we had a theoretically predicted L — T relation, cluster dimming would be accounted for by 
including a factor of -P 7 ^ 7 in the L — T relation, which would become: 

Similarly, the measured flux from a cluster in the 0.5-2 keV band, which was given by 
Eq. ( 2.24D , would become, 



L .5-2iW*) J<S(t+z) d f P T(f) 



\-'h^-„,r-)- fl 5 dfP T (f) 
However, since dimming mechanisms would mimic luminosity evolution, it is redundant to 



include the factor P 7 ^ 7 in Eq. (3.4) if the evolution specified by the parameter a is fit to 
observations. It is a remarkable fact that, although dimming does affect the appearance of 
clusters in these surveys, all these effects are absorbed into the z-dependence of the L — T 
relation. Thus, having measured this evolution, these surveys should determine the nature of 
dark energy independently of dimming. On the other hand, physics which would affect the 
formation of structure rather than the appearance of clusters, can be constrained with cluster 
counts without a theoretical understanding of the evolving L — T relation. 

Let us note here two important upshots of this fact: first, given a measured L — T 
evolution, implications of cluster counts are independent of dimming while supernovae clearly 
are not. Thus, as data sets expand, comparing these two will constrain any anomalous 
dimming of supernovae. Secondly, distant clusters tend to be brighter than they would have 
been in the absense of evolution. With a larger statistical sample of x-ray clusters and a better 
theoretical understanding of L — T evolution, this alone may be the strongest constraint on 
dimming mechanisms. 



- 13 - 



4. Results 



In this section we compare number counts found using the above model to the 400d ROSAT 
survey |2£(j. The 400d survey identified 242 optically verified X-Ray sources in the main 



survey. The search was done with a flux limit of 1.4 x 10~ 13 erg s~ x cm~ 2 and with a 
geometric sky coverage of 446.3 square degrees. In order to compare our theoretical number 
counts to ROSAT's data, we integrated redshift over a bin size of Az = 0.1. Since the 400d 
survey reported the error bars in their flux measurements, we estimated the error bars on 
galaxy cluster number counts by counting how many objects in a given redshift bin would lie 
below the flux limit when the measured flux is shifted downward by one standard deviation. 
In the cases where an X-ray source lied on the boundary of a redshift bin, it was counted 
in both of the adjacent redshift bins. To account for the scatter in the M — T and L — T 
relations, we included a log normal distribution in the effective L — M relation (L(M, z)): 



P L (lnL',z) 



: exp 



(ln(L')-ln(L(M,z))) s 
2<L 



(4.1) 



The effective selection probability P se /(M, z) of objects of mass M at redshift z is a convolution 
of the survey selection function P se l(f) with the distribution in luminosities inferred from the 
L-M relation P L {L, z): 



Psel(M,z) 



In L x (z) 



Pl(L',z)P Si 



JnL' 



I 4ird L (i 



din I/, 



where 



L x (z) = A7rd L (z) 2 f x 



(4.2) 



(4.3) 



is the lower limit on the luminosity at redshift z, corresponding to the flux limit f x of the 
survey, and the argument of P se \ is the flux expressed in terms of the luminosity and luminosity 
distance. We compared the best fit values of various observables as the assumed scatter, G\ n L, 
varied from 0.3 to 0.7 [35, 54]. The results are described below, and can also be seen in Fig. |l] 
and Fig. ||| 

The number of observed virialized objects in a redshift bin Az = 0.1 is then given by 
integrating the mass function over objects, weighted by P se i(M, z): 



N(> f x ,z,Az)=A_ 



geo 



z+Az 



z-Az JO 



dr dN ~ 

dzdMr(z) 2 - — P sel (M,z), (4.4) 



where A geo is the geometric sky coverage in steradians and f x is the flux limit of the survey, 
which feeds into P se i(M,z) as described above. 

The comoving volume element per steradian is, 



dV{z) = r{zf 



dr(z) 
dz 



dz, 



(4.5) 
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with 



and the Hubble parameter H(z) is given by Eq. (2.5). 

It is important to ensure that the lightest mass virialized object included in N given the 
assumed M — T and L — T relations is cluster sized and not smaller. Otherwise N contains 
smaller objects which are not included in the survey. This constrains the smallest z for which 
this formalism is valid. In our fits we only include clusters with redshift z > 2. The lightest 
mass virialized object that could have been observed by the 400d survey, with X-ray flux limit 
1.4xl0~ 13 erg/s/cm 2 , assuming Q, m = 0.3 and h = 0.73, is around M(0.2) = 1.6 x 1O 14 M , 
which is cluster size. 

4.1 Systematic Errors 

The only errors included in our fits are an estimate of the uncertainty in low flux cluster counts. 
There are in addition a number of theoretical uncertainties, some of which we studied by 
repeating our fits with different parameter choices. Larger normalizations of T15 in Eq. ( [2.17 ) 



would lead to a prediction of brighter clusters, and hence larger cluster counts. For a fixed 
data set, larger T15 would then translate into a measurement of less structure, corresponding 
for example to smaller Q m and/or erg. Similarly, a larger normalization for A§ in Eq. ( [2.18 ) 



would imply brighter clusters, with similar consequences to increasing T\§. To examine the 
uncertainty in predictions for cosmological and dimming parameters, we repeated our fits for 
typical determinations of T15 from simulations, differing by as much as 20%. The uncertainty 
in Aq is effectively equivalent to an additional uncertainty in T15 of around 5%. The power 
laws used in the L — T and M — T relation also have associated errors, and it would be useful 
to perform a more complete error analysis. Another source of error is our assumed power law 
spectrum used to calculate K-corrections, which is a worse approximation for low-temperature 
clusters than high-temperature clusters, but we expect that this is not a significant source of 
error. We have also checked that alternative redshift binning of the data does not significantly 
change our results. To examine the effect of scatter on our analysis, we reduced our assumed 
scatter from (J\ n L = 0.7 to 0.3 and found that the predicted number counts were reduced by 
nearly a factor of two, as can be seen below. This demonstrates the importance of correctly 
accounting for such statistical effects. 

4.2 Flux Limited Cluster Counts for Standard Cosmology 

Figure || shows our computed number counts and the 400d survey's observed number counts 
versus redshift. Three curves were drawn for different values of O m <jg, with T = 0.2 and 
T15 = 6 keV. 

Larger normalizations for the L — T relation (T15) lead to smaller predicted values of 
Q m and erg. Reducing the assumed level of scatter in the effective L — M relation leads to a 
decrease in the predicted number of dim objects observed. For a given set of observations, 
reducing the assumed scatter would then lead to a larger inferred amount of structure, i.e. 
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Figure 1: Number Counts versus redshift for different matter densities (f2 m ) and matter density 
fluctuation amplitudes (as), without photon-axion oscillations, with different levels of scatter in the 
L-M relation. The theoretical predictions correspond to T = 0.2 and T15 = 6 keV. (a) a\ n L = 0.3. (b) 
a\nL = 0.7. 



larger a s . For T 15 = 6 keV and T = 0.2, the best fit shifts from n m = 0.209 and a 8 = 0.923 
with a\ n i, = 0.3 to J7 m = 0.286 and erg = 0.731 with a\ n L = 0.7. 

Figure ^ shows our \ 2 analysis for different values of T15, T and a\ n L- We only include 
our estimated uncertainties in the 400d survey number counts in the statistics. Notice that 
for values of V between 0.1 and 0.2 and T15 between 5 and 6 keV, there is tension between our 
result and the best fit WMAP 3-year measurement. We are consistent with WMAP bounds 
if we assume large T15, small V and large (J\ n i- Our results are similar to earlier studies 
[||, [D], Hi, although Reiprich @ has found that the HIFLUGCS cluster survey is in still 
better agreement with the WMAP 3-year Q and COBE 4- year data [55|. Flux limited cluster 
counts can also be used to constrain other cosmological parameters, such as the equation of 
state parameter w (for example, M). 



4.3 Ineffectiveness of Flux Limited Cluster Counts for Dimming Mechanisms 

As we mentioned earlier, we cannot use cluster counts to constrain dimming mechanisms. This 
is not to say that dimming mechanisms do not affect cluster counts; indeed, dimming would 
lead to fewer clusters above the flux limit in any given survey. However, the effect of dimming 
would only be through a modification of the observed evolution of cluster luminosities, which is 
currently not well constrained theoretically. To gauge the effect that photon-axion oscillations 
could have on cluster counts we can assume some particular intrinsic L — T evolution, and 
examine the predicted number counts with and without oscillations. Figure ^ shows number 
counts versus redshift for different values of L^ ec , where Ld oc = 00 corresponds to no photon- 
axion oscillations and values of L^ec as low as 30 Mpc correspond to roughly 1/3 of the light 
lost. The curves correspond to O m = 0.3 and a 8 = 0.9 with T = 0.1 and T15 = 6 keV. We 
assumed here that the intrinsic luminosity evolution is specified by the parameters ( [2.22 ), 
although there is no theoretical justification for this. 
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h=0.73, r=0.2, o-| nL =0.3, and L dec ->oo h=0.73, r=0.2, a-| nL =0.7, and L dec ->oo 




0.2 0.22 0.24 0.26 0.28 0.3 0.2 0.22 0.24 0.26 0.28 0.3 



(c) (d) 

Figure 2: Confidence plot of Vl m and as for various choices of model parameters. The lines represent 
68%, 80%, 90%, and 95% confidence regions, (a) T = 0.1, cri„ L = 0.3. (b) V = 0.1, a lnL = 0.7. (c) 
r = 0.2, a XnL = 0.3. (d) r = 0.2, a lnL = 0.7 



The fact that the observed luminosity evolution is used as input in this analysis implies 



that the constraints on cosmological parameters from Section 4.2 are valid, independent of 
any dimming mechanism. As a consequence, cluster constraints on the equation of state 
parameter w can be compared with constraints from Type la supernova surveys, which would 
be affected by dimming mechanisms, as in Refs. Such a comparison would then provide a 
new test of the photon-axion oscillation model. It would also be interesting to compare with 



other constraints on photon-axion oscillations, for example from CMB spectral distortion [56]. 

4.4 Flux Limited Cluster Counts for Other Types of New Physics 

Although we do not attempt further analyses here, flux limited cluster counts are more 
suitable for constraining new physics that modifies structure formation as opposed to the 
apparent luminosity of clusters. There are many well-known examples of such possible new 
physics. These include possible new interactions in the dark sector and new light species 
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Figure 3: Number Counts versus redshift for Cl m — 0.3 and as = 0.9 including axion oscillations for 
a fixed intrinsic L-T relation, with different levels of scatter in the L-M relation, (a) cry n L = 0.3. (b) 
CTlnL = 0.7. 



that would wash out structure. Cluster surveys would also be useful in constraining such 
phenomena as late time phase transitions in the dark sector and other aspects of possible new 
gravitational dynamics. It would be straightforward to build these types of new physics into 
the formalism described above. 



5. Significance of Future Cluster Surveys 

It is important to recognize that studies of cluster evolution, while already interesting, will 
continue to develop. In this section, we briefly mention some future approaches that will 
enhance our knowledge of cluster growth, and note the impact of dimming. In particular, 
Sunyaev-Zeldovich (SZ) surveys such as the South Pole Telescope (SPT) |57j or the Atacama 
Cosmology Telescope (ACT) [gg] will establish catalogues of clusters which are unbiased in 
redshift - a crucial element difficult to achieve with X-ray surveys. The Dark Energy Survey 
(DES) H| will take advantage of the SPT survey, and include photometric redshifts, as well 
as lensing measurements of cluster masses, and other, independent tests of cosmology. The 
Large Synoptic Survey Telescope (LSST) |30| will provide masses of a huge set of clusters 
via weak lensing tomography. Future x-ray surveys can expand tremendously the statistics 
and knowledge of many of the uncertainties described in earlier sections, in particular the 



evolution of cluster properties [61]. 



One of the key difficulties in using X-ray surveys to extract cosmology, especially within 
the present context of photon- axion oscillations, is the indirect, and uncertain relationship 
of luminosity to temperature and temperature to mass. Future studies will mitigate these 
issues. 

SZ surveys will be a tremendous source of new information in the near future. For a 



review, see [62]. The SZ effect is a decrement AT in the CMBR given by the line of sight 
integral 

AT a T 



TcMBR rrieC 



dln e {l)k B T e {l) (5.1) 
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where cjt is the Thomson cross section, m e is the electron mass, n e is the electron number 
density and T e is the electron temperature. Clusters, with masses in excess of 1O 14 M , 
have sufficient gas densities and temperatures that large scale surveys are possible. The key 
feature of the SZ effect is the perturbation of the CMBR which is independent of redshift, 
and thus allows for a cluster sample, without concern of the selection issues associated with 
luminosity-weighted X-ray surveys described in earlier sections. 

The SZ effect is not proportional to mass alone, but to the electron pressure. Extracting 
the mass is a challenge, and of vital importance if these surveys are to provide precision 



limits on cosmology p3| , 64, |65|]. Techniques can involve self calibration |66|, 67, 64, |68| , 
measurements of the cosmic shear (as in the DES or LSST), or complementary measurements 
of the cluster x-ray temperature. 

As already noted, microwave studies, such as SZ surveys, should not be impacted by dim- 
ming mechanisms. Hence, the appearance and properties of clusters within such experiments 
should be a robust test of the dark energy properties. Similarly, surveys employing weak lens- 
ing will also determine mass and redshift properties will also be insensitive. Supernovae, on 
the other hand, acting as standard candles, are clearly impacted. As a consequence, studies 
of cluster growth are key tools of cosmology and measure a quantity distinct from that of 



supernovae, quite in the way envisioned by Wang, et al [21|. 
6. Conclusions 

We have reviewed some of the current models of structure formation and galaxy cluster 
dynamics relevant for comparing cluster surveys with models of particle physics and gravity. 
We compared predictions in the standard ACDM cosmology to the 400 Square Degree ROSAT 
galaxy cluster survey, and found that, with a relatively large assumed scatter in the relation 
between cluster mass and temperature, our analysis is consistent with the WMAP 3-year 
data. Earlier analysis of the HIFLUGCS cluster survey indicates even better agreement with 



CMB data |0|. 



We studied a model of cluster dimming by photon-axion oscillations, and found that 
a better theoretical understanding of cluster luminosity evolution is required before firm 
conclusions could be drawn regarding dimming mechanisms using cluster data. In particular, 
improvements in theoretical models and experimental measurements of the evolution of the 
luminosity-temperature relation may provide an important test of such mechanisms in the 
future. Moreover, we noted that the cosmological parameters extracted from cluster count 
surveys are independent of dimming mechanisms, given the measured L — T relation, in 
contrast to supernovae, and thus provide an independent test of such models. This suggests 
that such surveys should be folded into analyses such as in order to additionally constrain 
them. Although cluster counts are insensitive to dimming mechanisms, it is possible to 
impose relatively strong constraints on models of new physics that would affect structure 
formation as opposed to cluster appearances. This includes any physics that would alter 
the overall expansion rate of the universe, the existence of light species that would help to 
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wash out structure, and additional interactions in the dark sector that could either encourage 
or inhibit the growth of structure. The effects of new physics on structure formation can 
straightforwardly be built into the Press-Schechter formalism. 

The simple models that enter our fits rely on comparison to numerical simulations for 
justification, and deserve to be scrutinized. For example, a higher than typical normalization 
T 15 in the M — T relation, while not justified by simulations, could bring our fits still better in 
line with the WMAP 3-year data. Also, the level of scatter in the L — T and M — T relations 
deserves to be better analyzed, since a high assumed scatter brings our fits better in line with 
the WMAP fits. Despite the remaining uncertainties in this formalism, galaxy cluster surveys 
are increasingly ripe for their utilization in constraining new physics. Upcoming Sunyaev- 
Zel'dovich and weak-lensing cluster surveys hold the promise of a better determination of the 
cluster mass function, and should help to eliminate some of the current uncertainties in these 
techniques. 
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